function im_out=remap_res_im(im_in)


midpoint=ceil(size(im_in,2)/2);
[~,col_ind]=unique(round(sin(([1:size(im_in,2)]-midpoint)/midpoint*24/41.7*pi/2)*midpoint+midpoint));

im_out=im_in(:,col_ind);

